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ABSTRACT 


There  are  two  methods  which  have  been  developed  inde¬ 
pendently  for  computing  network  sensitivities.  Both  computa¬ 
tions  may  be  carried  out  in  the  frequency  or  in  the  time 
domains.  One  method  involves  the  analyses  of  two  networks  - 
the  original  and  its  mutually  reciprocal  adjoint.  The  second 
method  uses  a  sensitivity  model  for  the  circuit.  It  is  shown 
that  the  sensitivity  model  and  the  mutually  reciprocal  adjoint 
circuit  are  essentially  the  same;  the  sensitivity  model  being 
useful  for  calculating  single  parameter  sensitivity  in  the 
time  domain/  the  adjoint  circuit  being  useful  for  calculating 
sensitivity  for  several  parameters  in  the  frequency  domain. 
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I .  INTRODUCTION 


The  subject  of  sensitivity  is  one  of  the  oldest  areas 
where  electrical  engineers  have  made  extensive  studies.  A 
great  amount  of  research  work  has  been  done  and  in  almost 
every  technical  journal  articles  have  been  published  concerning 
different  approaches  to  network  sensitivities.  This  subject 
received  new  impetus  when  the  digital  computer  was  made 
available  to  almost  every  electrical  engineer.  With  the  aid 
of  the  computer,  different  approaches  to  computer  aided 
circuit  design  have  been  outlined  using  -  sensitivity  models  or 
the  mutually  reciprocal  adjoint  network. 

Starting  with  th'e  definition  of  "Interreciprocity,"  S.  W. 
Director  and  R.  A.  Rohrer  [Ref.  1-4]  developed  the  idea  of 
automated  network  design  and  sensitivity  calculations  for 
linear,  time  invariant,  and  later  for  nonlinear,  time  variant 
circuits,  using  the  reciprocal  adjoint  network.  The  calcu¬ 
lations  can  be  carried  out  in  the  frequency  or  in  the  time 
domains,  although  time  domain  calculations  are  involved. 
Computations  of  the  sensitivity  due  to  changes  in  all  network 
parameters  require  the  simultaneous  analysis  of  two  networks, 
which  is  easily  accomplished  with  the  aid  of  a  digital  computer. 

The  other  approach  to  network  sensitivity  makes  use  of 
sensitivity  models  as  developed  and  published  by  J.  V.  Leeds 
and  G.  I.  Urgon  [Ref.  5] .  These  results  were  extended  later 
by  S.  R.  Parker  [Ref.  6]  to  nonlinear  time-variant  circuits. 
Using  sensitivity  models  the  changes  of  an  output  quantity 


due  to  variations  of  one  circuit  parameter  are  easily  achieved, 
for  complicated  networks  the  computations  of  the  sensitivity 
due  to  changes  in  all  network  parameters  are  more  involved. 

The  computations  are  carried  out  in  the  frequency  and  the  time 
domains  equally  well. 

It  is  the  subject  of  this  thesis  to  show  that  both 
approaches  to  network  sensitivity  are  not  independent.  First 
a  careful  review  of  the  mutually  reciprocal  adjoint  network 
is  given.  As  a  new  result  a  topological  relationship  between 
the  original  network  and  its  adjoint,  including  dependent 
sources  and  independent  sources,  is  presented  as  noted.  After 
that  the  relations  and  transitions  between  the  adjoint  network 
approach  to  sensitivities  and  the  sensitivity  model  are  shown. 
Finally  the  advantages  and  disadvantages  of  both  methods  are 
discussed. 


II.  THE  MUTUAL  INTERRECIPROCAL 
ADJOINT  NETWORK 


The  interreciprocity  property  of  an  original  network  N 
and  its  mutual  adjoint  network  N  is  an  important  extension 
of  the  reciprocity  theorem  used  for  computation  of  multi¬ 
parameter  sensitivities  and  automated  network  design.  As 
defined  by  Director  and  Rohrer  [Ref.  1]  the  properties  are 
summer ized  in  the  following  paragraph. 

A.  THE  ADJOINT  NETWORK 

For  any  general  network  N  containing  arbitrary  multi¬ 
terminal  or  two-port  -elements  with  parametric  representation 

»*» 

(lumped  parameters) ,  there  exists  an  adjoint  network  N  which 
has  the  same  topology,  but  not  necessarily  the  same  element 
types,  in  corresponding  branches. 

1.  The  Linear  Time  Invariant  Case 

Director  and  Rohrer  [Ref.  2  and  3]  developed  the 
adjoint  network  N  as  being  identical  to  the  original  network 
with  the  following  exceptions: 

a)  All  gyrators  in  N  with  gyration  ration,  a,  become  gyrators 
in  N  with  gyration  ration,  -a,  (polarity  reversed) . 

b)  All  voltage  controlled  voltage  sources  in  N  become  current 
controlled  current  sources  in  N  and  voltage  amplification  factor, 
y,  becoming  current  amplification  factor,  -p. 

c)  All  current  controlled  current  sources  in  N  become  voltage 

•v 

controlled  voltage  sources  in  N  with  controlling  and  controlled 
branches  reversed  in  N  and  current  amplification  factor,  h, 
becoming  voltage  amplification  factor,  -h. 
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d)  Ail  voltage  controlled  current  sources  and  current 
controlled  voltage  sources  have  their  controlling  and  controlled 
branches  reversed  in  H. 

e)  All  independent  sources  are  set  to  zero.  For  computations 
of  sensitivities  of  network  functions ,  excitations  with  unity 
sources  at  specified  terminals  are  explained  in  a  later 
paragraph . 

f)  In  the  frequency  domain  no  changes  occur  in  the  excitation 
of  the  two  networks.  In  the  time  domain,  time  in  the  adjoint 
network  runs  backwards. 

If  the  two-port  coupling  elements  such  as  transformers, 
gyrators,  and  dependent  sources  are  defined  by  algebraic 
relations  among  their  port  voltages  and  currents,  then  these 
relations  can  be  summerized  as  shown  in  Fig.  1.  The  ideal 
transformer,  the  voltage  controlled  voltage  source,  and  the 
current  controlled  current  source  are  described  by  the  hybrid 
matrix.  The  gyrator  is  expressed  either  by  the  impedance  or  the 
admittance  matrix.  The  current  controlled  voltage  source  is 
defined  by  the  impedance  matrix,  and  the  voltage  controlled 
current  source  by  the  admittance  matrix.  In  Fig.  1  the  first 
subscript  is  defined  as  follows: 
i  -  input  branch 
o  -  output  branch 

The  second  subscript  denotes  the  kind  of  two-port  element  and 
is  defined  as  follows: 

y  -  voltage  dependent  voltage  source 
h  -  current  dependent  current  source 


/ 
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g  -  voltage  dependent  current  source 
r  -  current  dependent  voltage  source 
n  -  ideal  transformer 
c  -  gyrator 

The  second  subscript  is  omitted  in  the  matrix  representation 
but  will  be  used  later. 

Let  the  voltages  and  currents  belonging  to  branches 
in  the  mutual  reciprocal  adjoint  network  be  defined  by 


and  1 


xy 


respectively.  The  subscripts,  xy,  are  explained  as  used  later 
on. 

The  transformation  of  all  passive  circuit  elements 
from  the  original  network  N  into  its  corresponding  adjoint  N 
can  then  be  summerized  as  shown  in  Fig.  2a  and  b.  These 
transformations  are  valid  for  any  linear  time  invariant  network. 
As  stated  before,  for  sensitivity  calculations  all  independv..-* . 
sources  in  the  original  and  its  adjoint  network  are  set  to 
zero. 

2.  The  Linear  Time  Variant  Case 

The  adjoint  network,  N,  of  the  original  network  N 
for  the  linear  time  variant  case  is  defined  by  Director  and 
Rohrer  [Ref.  1]  as  follows: 

a)  All  time  invariant  elements  of  N  become  elements  in  N 
as  described  in  the  previous  paragraph. 

b)  All  time  varying  resistors,  gyrators,  transformers,  and 
controlled  sources  of  N  are  time  varying  in  N.  The 
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transformations  are  according  to  the  rules  governing  the 
corresponding  time  invariant  elements. 

c)  Time  varying  capacitors,  C(t),  of  N  become  time  varying 
capacitors,  C(x),  of  N  shunted  by  a  time  varying  conductance, 

G (t) ,  in  mhos  equal  to  the  value  of  the  time  derivative  of 
the  capacitor. 

d)  Time  varying  inductors,  L(t),  of  N  become  time  varying 
inductors,  L(t),  of  N  in  series  with  a  time  varying  resistance, 
R(t) ,  in  ohms  equal  to  the  derivative  of  L(t)  with  respect  to 
time. 

e)  Time  varying  coupled  inductors  and  their  adjoint 
equivalent  are  shown  in  Fig.  3c. 

f)  In  the  time  domain  calculations,  time  in  N  runs  backwards 
relative  to  time  in  N.  If  the  initial  time  is  defined  by  tQ, 
the  final  time  by  t^,  and  the  running  time  in  N  by  t,  then  the 
time  in  N  is  given  by 

X  -  t  +  t-  -  t  ;  t  <  t  <  t_  (1.1) 

of  o  f 

g)  The  adjoint  network  for  frequency  calculations  is  identical 
to  the  adjoint  network  in  the  time  domain,  except  there  is  no 
backward  running  time.  For  the  sensitivity  calculations  the 
network  analyses  of  both  networks  has  to  be  carried  out  at 

each  frequency  point  simulataneously . 

h)  All  independent  sources  are  set  to  zero.  The  network 
excitations,  for  specific  sensitivities  of  a  network  function, 
are  discussed  later. 
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Figure  1:  Two-Port  Elements  and  their  Matrix 
Characterization 

a)  Voltage  Dependent  Current ' Source 

b)  Current  Dependent  Voltage  Source 

c)  Voltage  Dependent  Voltage  Source 

d)  Current  Dependent  Current  Source 

e)  Ideal  Transformer 

f)  Gyrator 
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B,  THE  INTERRECIPROCAL  THEOREM 

The  interreciprocal  property  of  an  original  network  N 
and  its  adjoint  N  as  defined  by  Director  and  Rohrer  [Ref.  1] 
is  a  generalisation  of  the  reciprocity  theorem.  This  extension 
applies  to  a  network  and  its  adjoint  consisting  of  resistors,, 
capacitors,  inductors,  coupled  inductors,  transformers,  gryators 
and  controlled  sources.  The  reciprocity  theorem  defines  a 
network  to  be  reciprocal  if  it  has  the  following  property: 

If  an  excitation  E  is  applied  at  one  pair  of  terminals 
in  N  and  a  response  I 2  is  measured  at  some  second  pair  of 
terminals  of  the  same  network,  interchanging  the  points  of 
excitation  and  response,  keeping  E  the  same,  does  not  change 
the  response  I2  at  the  original  port  (Fig.  4) . 

An  original  n-port  network  N  and  its  n-port  adjoint  N 
are  said  to  be  interreciprocal  if  the  following  conditions  are 
satisfied: 

Considering  first  the  frequency  domain  case.  For  any 

excitation  E^(s)  at  some  terminal  pair  k  of  the  original 

network  N  the  response  at  another  terminal  pair  n  is  In(s) . 

The  excitation  at  all  other  ports  is  zero.  Exciting  the  adjoint 

network  at  terminal  pair  n  with  the  sour  V  (s)  such  that 

n 

Vn(s)  =  Er(s)  (1.2) 

yields  the  response 

Ik(s)  =  In(s)  (1.3) 

at  terminal  pair  k  of  N. 
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For  time  domain  considerations  let  the  excitation  voltage 
at  port  k  of  N  be  e^t)  which  forces  a  current  response  in(t) 
at  the  terminal  pair  n.  Exciting  the  adjoint  at  its  port  n 
with  a  voltage  source 

Vn(t)!  “  ek(t)  '*  T  =  t0  +  tf-t  (1.4) 

t=t 

yields  the  current  response  i^Ct)  at  terminal  pair  k  of  N 
such  that 

±k(t)|  =  1  (t)  (1.5) 

K  t=T  n 

The  voltage  and  current  excitations  at  all  other  terminal  pairs  are 
zero  (Fig.  5).  If  these  conditions  (1.2  through  1.5)  apply 
to  all  possible  pairs  of  terminals  of  both  networks  then  they 
are  said  to  be  interreciprocal. 

One  sufficient  condition  for  an  original  network  and  its 
adjoint  to  be  interreciprocal  is  that  the  circuit  consist  of 
linear  time  invariant  parameters  only. 

C.  DEFINITION  OF  THE  ADJOINT  NETWORK  IN  TERMS  OF  TOPOLOGICAL 

RELATIONSHIPS 

The  interreciprocity  theorem  applied  to  a  network  and  its 
adjoint,  implies  certain  restrictions  on  the  transformation  of 
elements  from  one  circuit  to  the  other.  It  results  in  very 
strict  relationships  between  the  original  network  and  its 
adjoint.  As  discussed  by  Parker  and  Barmes  [Ref.  7]  the  branch 
relations  of  the  original  network  can  be  expressed  by  the 
following  matrix  equation 
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(2.1) 
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where  and  are  defined  as  link  voltages  and  currents, 
respectively,  v2  and  i2  are  branch  voltages  and  currents, 
respectively.  They  are  defined  by  the  following  vectors 


VS 

XS 

< 
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(2.2) 


(2.3) 


where  the  subscripts  denote  the  following: 

C  ~  tree  branch  capacitances 
G  -  tree  branch  conductances 
T  -  tree  branch  (excess)  inverse  inductors 
S  -  link  (excess)  susceptances 
R  -  link  resistors 
L  -  link  inductors 

Independent  sources  consist  of  voltage  and  current  sources  as 
defined  by  the  following  vectors: 


V 

e  = 

SR 

and  j  = 

^G 

eL 

(2.4) 
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where  e^,,  eR,  and  are  independent  voltage  sources  contained 
in  fundamental  loops  defined  by  susceptances,  resistors/  or 
inductors,  respectively.  jc,  jG,  and  jp  are  independent  current 
sources  associated  with  fundamental  cutsets  defined  by  capacitors, 
conductors  or  inverse  inductors,  respectively.  The  latter 
can  be  neglected  in  the  adjoint  network  because  all  independent 
sources  are  set  to  zero.  Only  indpendent  sources  relevant  to 
the  computation  of  the  sensitivity  are  inserted.  These  appear 
as  unity  sources  in  a  reciprocal  manner  in  the  original  network 
and  its  adjoint  as  discussed  later. 

The  two  partitioned  square  matrices  (2.1)  contain  various 
elements  as  follows: 


Mx  =  f (u1  +  I) ] 
Y1  =  [gl] 


^1  ^rl^ 


(2.5a,b,c,d) 


J;L  =  [(hx  +  I)] 


and 


Z  =  [r]  M  =  [(M  +  F)] 

L  =  [ (h  -  Ffc) ]  Y  =  [g]  (2.6a,b,c,d) 


where 

consists  of  voltage  dependent  voltage  sources  plus 
the  identity  matrix 

Z^  contains  current  dependent  voltage  sources 

is  composed  of  voltage  controlled  current  sources 
consists  of  current  dependent  current  sources  plus 
the  identity  matrix 

Z  contains  the  remaining  current  dependent  voltage  sources 
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M  contains  the  other  part  of  voltage  dependent  voltage 
sources  plus  partitioned  topological  F  matrix  elements 
L  consists  of  further  current  dependent  current  sources 
plus  the  negative  transposed  elements  of  F 
Y  is  composed  of  all  remaining  voltage  controlled 

current  sources.  Using  (2.2)  (2.3)  (2.5),  and  (2.6) 
in  (2.1)  with  e  =  j  =  0  yields 
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where  the  double  subscripts  indicate  the  kind  of  elements  between 
which  the  dependency  exists. 

The  basic  transformation  from  the  original  to  the  adjoint 
contains  no  changes  for  all  passive  circuit  parameters.  Voltage 
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controlled  voltage  sources  become  current  controlled  current 
sources  with  current  amplification  factor ,  -y,  with  the  roles 
of  controlling  and  dependent  branches  reversed.  This  operation 
corresponds  to  the  following  matrix  manipulation 


•  1 

r  i  i 
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i  i 
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t 

— 7 - 
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i  t 

1 

L  hi  j 

and 


(2.8) 


(2.9) 


A  similar  transformation  holds  for  current  dependent  current 
sources  so  that 


-h 


1~ 

h 

l 


(2.10) 


and 


t 


(2.11) 


Voltage  dependent  current  sources  and  current  dependent  voltage 
sources  with  amplification  factor,  g,  and  amplification  factor, 
r,  respectively  remain,  but  in  both  cases  the  roles  of  depending 
and  controlling  branches  are  reversed.  This  operation 
corresponds  to  the  transposition  of  the  corresponding  submatrices 
so  that 
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(2.12) 
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Equation  (2,1)  in  partitioned  form  and  omitted  independent 
sources  gives 
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Using  the  matrix  transformation  as  shown  in  (2.8)  through  (2.13) 
in  (2.14)  yields  the  branch  relations  of  the  mutual  reciprocal 
adjoint  network 
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Rearranging  and  using  (2.5)  and  (2.6)  in  (2.15)  gives  the  final 
result 


(2.16) 
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Compared  with  (2.1)  shows  an  easy  and  compact  relationship 
between  the  original  network  and  its  adjoint. 


D.  USE  OF  THE  ADJOINT  NETWORK  FOR  SENSITIVITY  CALCULATIONS 

R.  A.  Rohrer  and  S.  W.  Director  [Ref.  1,  2,  and  4]  have 
shown  that  the  sensitivity  of  any  network  function  with  respect 
to  changes  of  one  or  all  network  parameters  requires  the 
analysis  of  the  original  circuit  and  its  mutual  reciprocal 
adjoint. 

If  any  network  function  is  denoted  by  H(jw),  then  inserting 
current  or  voltage  sources  of  one  ampere  or  one  volt,  respec¬ 
tively,  at  particular  ports  of  N  and  N,  excited  in  a  reciprocal 
manner,  H(jw)  becomes,  a  network  function  as  shown  in  Table  1. 


TABLE  1.  Definition  of  the  Network  Emmittance  Function. 


- - - ■ — ’ — 

u 

n  ijw) 

Terminal  Conditions 
Original  Network 

N 

Terminal  Conditions 
Adjoint  Network 

M 

port  k 

port  1 

port  k 

port  1 

Driving  Point  Impedance 
at  port  k 

Current 
source 
(1  amp) 

Open  or 
Short 

Current 
source 
(1  amp) 

Open  or 
Short 

Driving  Point  Admit¬ 
tance  at  port  k 

Voltage 
source 
(1  volt) 

Open  or 
Short 

Voltage 
source 
(1  volt) 

Open  or 
Short 

Transfer  Impedance 
between  port  k  and  1 

Open 

Current 
source 
(1  amp) 

Current 

source 

(lamp) 

Open 

Transfer  Admittance 
between  port  k  and  1 

Short 

Voltage 
source 
(1  volt) 

Voltage 
source 
(1  volt) 

Short 

Current  Transfer  Ratio 

Current 
source 
(1  amp) 

Short 

Short 

Voltage 
source 
(1  volt) 

Voltage  Transfer  Ratio 
between  port  k  and  1 

Voltage 
source 
(1  volt) 

Open 

Open 

Current 
Source 
(1  amp) 
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The  normalized  sensitivity  (due  to  the  insertion  of  unity 
current  or  voltage  sources)  of  any  network  function  (as  defined 
in  TABLE  1)  with  respect  to  all  element  types,  is  obtained  in 
terms  of  voltage  and/or  current  responses  in  the  corresponding 
branches  of  N  and  N.  The  sensitivities  are  defined  in  TABLE  2. 


TABLE  2,  Sensitivities  of  a  Network  Function 


Variable  Network 
Parameter 

Sensitivity  of 

Network  Function 

Resistances 

||  -  —  IR<jw)-IR(jw) 

Conductances 

If  =  VG(jw)-VG(jw) 

Inductances 

Reciprocal  Inductances 

~  =  1/ ( jwVp  (jw)  *Vp(jw)  ] 

Capacitances 

Elastances 

||  =  — 1/ (jwlg (jw) • Ig ( jw) ] 

Transformers  (turns  ratio 
n:l) 

IS  -  »on(3w)'Vin(5w>+Vi„(ju>- 

xon ( jw) > 

Gyrators  (gyration  ratio: 
a) 

|f  -  ^Uia(jw)-Iocl<jw)-Ioa(jw>- 

v»> 

Voltage  controlled  voltage 
sources  (voltage  ampli¬ 
fication  ratio:  ]i) 

If  =  -viu(jw)-i0M(jw) 

Voltage  controlled  current 
sources  (transconductance: 

g) 

If  =  Vig(jw)-Vog(jw) 

Current  controlled  current 
sources  (current  ampli¬ 
fication  ratio:  h) 

IS=  Iih(iw),voh(» 

Current  controlled  voltage 
source  (transresistance:  r) 

If  =  -Iir(jw)-Ior(jw) 

III.  SENSITIVITY  MODELS 


The  sensitivity  model  approach  to  network  sensitivities 
for  linear  circuits  was  developed  by  J.  V.  Leeds  and  G.  I. 

Urgon  [Ref.  5]  based  upon  an  idea  first  presented  by  R. 

Tomovic  [Ref.  8] .  These  results  were  extended  to  nonlinear 
circuits  by  S.  R.  Parker  [Ref.  6].  In  general/  the  sensitivity 
model  is  topologically  identical  to  the  original  circuit.  All 
independent  sources  are  reduced  to  zero.  An  excitation  voltage 
or  current  source,  depending  on  the  variable  parameter,  x,  has 
to  be  placed  in  series  or  in  parallel  with  x,  in  such  a  direction 
as  to  oppose  the  normal  current  flow  in  that  branch.  The  value 
of  that  source  depends  upon  the  current  or  voltage  response 
of  the  branch  of  x  in  the  original  network.  The  responses  of 
the  sensitivity  model  are  in  turn  the  required  sensitivity 
function. 

A.  DEFINITION  FOR  THE  LINEAR  CASE 

For  the  different  element  types  the  sensitivity  model 
equivalent  element  and  its  corresponding  excitation  is 
summarized  in  Table  3  as  taken  from  S.  R.  Parker  [Ref.  6]. 
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TABLE  3.  Sensitivity  Models  and  their  Excitations  for  the 
Linear  Network 


c)  CAPACITIVE  ELEMENT 


d)  INDUCTIVE  ELEMENT  (V7ITH  MUTUAL  COUPLING) 


3a 


lleq 

"  L11 

Vlleq  “ 

diL 

dt 

12eq 

-  L12 

V12eq 

dij, 

dt 

21eq 

=  L21 

V21eq  = 

diL 

dt~ 

22eq 

=  L22 

v  = 

22eq 

dip 

dt 

VL 

d 

dt 

^lii1l  *  Li21r) 

vr 

d 

dt 

(L21XL  +  L221r) 

* 


B.  DEFINITION  FOR  THE  NONLINEAR  CASE 

For  the  nonlinear  circuit  the  sensitivity  model  equivalents 
and  their  excitations  are  summerized  in  TABLE  4  as  taken  from 
S.  R.  Parker  [Ref.  6]  . 


TABLE  4,  Sensitivity  Model  and  Their  Excitation  for  the 
Nonlinear  Case 


a)  RESISTIVE  ELEMENT 
3V. 


R_ 


3a 


•-wvw — 

3i„  R  v — A 


3a 


■R  eq 


R 


'eq 


3vR(iR,a) 


eq 


IT 


'eq 


R 


3Vva> 

3a 


b)  CONDUCTIVE  ELEMENT 
3V_ 


3a 


VR  “  vR(Va) 


iG  =  iG(vG,a) 


C)  CAPACITIVE  ELEMENT 

3vc 


3VC 

'3a 


* - ►- 


_^C 

3a 

3QCvc,a) 

3vl 


eq  ( 


'eq 


'eq 


3 

3t 


3Q(vr,a) 

( - b - ) 

1  3a  1 


eq 


9Q  (vp/a) 

r  t - b - ) 

eq  3a  1 


LC  at  Q(vC'a) 
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d)  INDUCTIVE  ELEMENT  (WITH  MUTUAL  COUPLING) 


VL  3t  ^ll(lL'a)  +  ^12  (xr^a)  -1 


vr  3t  ^21^L/C1^  +  ^22(ir,a)  J 
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IV.  RELATIONSHIP  BETWEEN  THE  SENSITIVITY  MODEL 
AND  THE  MUTUAL  ADJOINT  NETWORK  APPROACH  TO  SENSITIVITY 

To  compare  the  two  methods  it  is  first  shown  how  network 
sensitivities  are  obtained  using  Tellegen's  theorem  in  conjunc¬ 
tion  with  an  original  network,  the  mutual  reciprocal  adjoint, 
and  the  augmented  original  network.  Following  this  derivation, 
sensitivity  models  are  shown  to  be  a  special  case  of  the 
mutual  adjoint  network. 

A.  PROOF  OF  NETWORK  SENSITIVITIES  USING  THE  ADJOINT  NETWORK 
AND  TELLEGEN'S  THEOREM 

In  chapter  II  it  was  stated  that  the  sensitivity  of  a 

network  function  is  obtained  using  the  response  of  an  original 

network  and  its  adjoint,  A  proof  is  presented  now. 

Consider  the  network  of  Fig.  6a,  excited  with  a  voltage 

source  E  at  port  1.  At  port  2  the  voltage  response  is  V„. 

g  ^ 

Fig.  6b  represents  the  same  circuit  with  all  of  the  elements 
augmented.  It  is  excited  with  an  identical  voltage  source  E^ 
at  port  1.  The  voltage  response  at  port  2  is  Vj  +  ^V2* 

6c  represents  the  adjoint  of  the  original  circuit  excited  in 
reciprocal  manner.  In  Fig.  6,  Xa,  represents  any  kind  of  one- 
port  passive  network  parameter.  To  apply  Tellegen's  theorem, 
the  port  voltages  and  currents  of  the  augmented  original  circuit 
and  the  adjoint  network  are  tabulated  as  follows: 
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Figure  6:  Reciprocal  Two  Port  Networks 

a)  Original  Network ,N 

b)  Augmented  Original  Network 

c)  Mutual  Reciprocal  Adjoint  Network, N 


Aug.  Cria.  Network 

Rec,  Adioint  Network 

Voltaqe  Current 

Voltaqe 

current. 

Port 

Is 

E  I, 

9  1 

+  AX^ 

0 

h  =  ™v2 

Port 

2: 

V2  +  AV2 

0 

V2 

I-  =  E 

2  g 

int . 

netw. : 

v«  +  4Va  *« 

+  AI 

a 

V 

a 

Xa 

Multiplying  and  adding  the  corresponding  terms  as  shown  above 
yields 

-EgV2  +  <Va  +  4V*a  +  EgV2  +  iV2Eg  =  0  (4-1) 

(Io,  +  *V;.  "  0  (4-2) 


Equating  (4.1)  and  (4.2)  and  rearranging,  results  in  the  basic 
expression  from  where  the  proof  starts  for  different  kinds  of 
network  parameter. 


(I  +  A I  )V 
a  a  a 


-  (V  +  AV  )I  -  E  V_  =  0 
a  a  a  g  2 


(4.3) 


1 .  Passive  Network  Parameters 

The  proof  is  presented  for  impedance,  inductive,  eind 
capacitive  parameters  only. 

The  constraints  for  the  impedance  case  are 


V  =  Z  I 
a  a  a 


V  -  Z  I 
a  a  a 


AV  =  AI  Z  +  I  AZ 
a  a  a  a  a 


(4.4) 
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Substituting  (4.4)  in  (4.3)  and  rearranging  gives 


iV2E9 


-I  I  AZ 


a  a  a 


(4.5) 


Letting  the  excitation  be  a  unit  voltage  or  current  source, 
respectively,  £  =  1,  leads  to 


AV-  =  -I  I  AZ 
2  a  a 


a 


(4.6) 


or 


-1  I 


a  a 


(4.7) 


Equation  (4.7)  gives  the  sensitivity  of  the  output  voltage 
with  respect  to  changes  in  one  impedance  parameter.  Multiplica¬ 
tion  of  the  current  through  the  variable  impedance  in  the 
original  network,  Ia,  and  the  current  through  the  corresponding 
parameter  in  the  adjoint  network,  1^,  is  done  conveniently  in 
the  frequency  domain. 

For  the  capacitive  parameter  the  constraints  are 


dv. 


a 


=  C 


a 


a 


dt 


=  jwC  V 
J  a  a 


dVa  . 

I  =  C  - — -  =  jwC  V 
a  a,.  J  a  a 
dt 


(4.8) 


AI  =  jwAV  C  +  jwV  AC 
a  J  a  a  J  a  a 


Substituting  these  constraints  into  (4.3)  and  solving  for 


AV0E  results  in 
2  g 

AV0E  =  jwV  V  AC 
2  g  J  a  a  a 


(4.9) 
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Letting  the  excitation  source  E^  equal  to  one  gives 


Tre  -  jwV  V 
AC^  J  a  a 


(4.10) 


Equation  (4.10)  is  the  sensitivity  of  the  output  voltage  of  the 
original  network/  N/  with  respect  to  perturbations  in  one 
capacitive  network  parameter. 

Finally  the  derivations  for  changes  in  an  inductive 
element  are  shown.  The  auxiliary  equations  are 


V  =  L  — -  =  jw  L  I 
a  adt  J  a  a 


V  =  I,  — -  =  jw  L  I 
a  <xdt  a  a 


(4.11) 


AV  =  jwl  AL  +  jwAI  L 
a  J  a  a  J  a  a 


Substituting  (4.11)  into  (4.3)  and  solving  for  AV„E  gives 

/  g 


AV0E  =  -jwl  I  AL 
2  g  J  a  a  a 


(4.12) 


Assuming  the  excitation  sources,  E^,  equal  to  one  gives 


=  -jwl  I 
AL^  J  a  a 


(4.13) 


Equation  (4.13)  gives  the  incremental  changes  in  output  voltage 
due  to  variations  in  one  inductive  element  in  the  original 
network. 

These  results  agree  with  the  given  relations  in 
TABLE  2,  developed  by  Director  and  Rohrer.  For  better  comparison 
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the  corresponding  expressions  to  equation  (4.7)/  (4,10)/  and 
(4.13)  are  repeated  here 

3R  =  ~IR^w)IR(jw) 

§§  =  jwVc(jw)Vc(jw) 

It  =  -jwiL(jw)iL(jw) 

If  the  sensitivity  of  the  output  voltage  depends  on 
variations  of  all  network  parameters,  the  increments  are  added, 
applying  the  principle  of  superposition.  The  summations  are 
taken  over  all  corresponding  network  parameters.  (4.7),  (4.10) 
and  (4.13)  then  become 

AV~  =  "El  I  AZ 
2  a  a  a 

AV~  =  EjwV  V  AC  (4.14a,b,c) 

2  aJ  a  a  a 

AV~  =  -Ejwl  I  AL 
2  J  a  a  a 

a 

2 .  Dependent  Sources 

As  an  example  for  all  four  kinds  of  dependent  sources, 
the  derivation  for  the  voltage  dependent  voltage  source  is 
presented.  The  proof  for  the  three  others  is  quite  similar.  In 
Fig.  7a  the  original  network,  excited  by  a  voltage  source  E 

y 

at  port  1  and  its  adjoint  (Fig.  7b),  excited  in  a  reciprocal 
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Port  1  port  2 


portal  port  2 


port  1  port  2 


Figure  7:  Reciprocal  Two  Port  Networks  Containing 
Voltage  Dependent  Voltage  Sources 

a)  Original  Network, N 

b)  Mutual  Reciprocal  Adjoint  Network, N 

c)  Augmented  Original  Network 


41 


manner ,  are  shown.  Fig.  7c  represents  the  augmented  original 
circuit  excited  by  the  same  voltage  source,  E  .  In  Fig.  7c 
the  incremental  voltage  change,  AVq^,  is  defined  as  follows 


AV  =  yAV.  +  V.  Ay 
oy  H  i  y  ry 


{4.15) 


The  ports  and  the  internal  voltages  and  currents  for  the 
augmented  original  and  the  adjoint  network  are  tabulate  and 
then  Tellegen's  theorem  is  applied 


Aug .  Orig 

.  Network 

Rec.  Adjoint  Network 

Voltage 

Current 

Voltage 

Current 

Port  1: 

E 

g 

h  +  AIi 

0 

.  h  “  ‘v2 

Port  2: 

Viv2 

0 

*2  "  Eg 

Internal 

Controlled 

** 

Side  : 

(Vi„+4Vo  > 

0 

viy 

tiIoy 

Dependent 
Side  : 

(I  +AI  ) 
oy  oy 

0 

Ioy 

Multiplying  the  inner  and  the  outer  columns  as  shown  above 
yields 

Outer  Product: 


"Va  +  (V2  ♦  AV2)Eg  +  J  <vi(1  ♦ 


+  Z  (V  +  AV  +  AV  )I  =  0 
iy/oy  oy  OU  oy  °y 


(4.16) 


Inner  Product: 


I  =0  (4.17) 

iy/oy 
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V0  =  H*E 
2  g 

then  for  the  augmented  network 
(V2  +  AV2)  =  (H  +  AH)  Eg 
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and 


AV0  =  AH*E  (4.22) 

2  g 

Substituting  (4.22)  into  (4.21)  yields 

AH*E2  =  E  “V.  I  Ay  (4.23) 

9  iy/oy  ly  oy 

Since  H  is  a  function  of  p,  an  incremental  change  in  the  voltage 
transfer  function  with  respect  to  the  voltage  amplification 
factor,  y,  is  given  by 


AH  =  E  {—)  Ay  (4.24) 

iy/oy  y 

Comparing  equation  (4.24)  with  the  rearranged  equation  (4.23) 
yields 


V  I 

3H  _  iy  oy 

ay 


(4.25) 


Letting  the  excitation  voltage  and  current  source,  E^,  equal 
to  one,  gives  the  voltage  transfer  function  sensitivity  with 
respect  to  the  voltage  amplification  factor,  y,  as  the  product 
of  the  controlling  branch  voltage  in  the  original  circuit  and 
the  dependent  branch  current  in  the  mutual  reciprocal  adjoint 
network.  This  proves  the  stated  result  of  TABLE  2  which  is 
repeated  here  for  convenience 


DH 

9y 


(jw)  I 


oy 


(jw) 
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B.  TRANSITION  BETWEEN  APPROACHES 

The  derivation  of  the  sensitivities  of  a  voltage  transfer 
function  in  the  previous  paragraph  was  carried  out  in  the 
frequency  domain.  The  sensitivity  model,  as  stated  in  Chapter 
III,  is  given  in  the  time  domain.  To  use  the  derived  equations 
in  the  time  domain  requires  further  interpretation.  For  the 
original  circuit  (Fig.  6a)  the  voltage  transfer  function  was 
defined  as  follows 


H  (s) 


V2(s) 

Vs> 


The  sensitivity  of  the  output  voltage  due  to  changes  of  any 
kind  of  passive  network  parameters,  u  ,  is  then  given  by 

•  ~  Ct 


3H(s,u  ) 

iV2<s’V  -  ys) Auc 


{4 .26) 


3u 


a 


where  the  parameter,  u^,  is  itself  a  function  of  s  and  x 
such  that 


ua  =  ua(s'x)  (4.27) 


Substituting  (4.27)  into  (4.26)  and  applying  the  chain  rule 
yields 


3H(s,u  )  3u 

AV,(s,u  )  =  E  (s )  •  ( - £_)(— “)Ax  (4.28) 

^  a  g  3u  3x 

a 


Using  for  impedance  type  parameters  equation  (4.7)  in  conjunction 
with  (4.22)  and  substituting  into  (4.28)  leads  to 


II  3u 

AV~  (s,x)  =  -2-2L  (— )Ax  (4.29) 

E  (s)  3x 
9 
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Making  the  assumption  that  E  (t)  is  a  unit  impulse,  u  =  Z, 

g  ct 

and  x  -  Z,  the  transition  into  the  time  domain  of  AV2  (s ,x) 
is  given  by  the  convolution  of  the  resistor  current  in  the 
original  network  and  the  current  through  the  corresponding 
resistor  in  the  adjoint  network.  Then  the  sensitivity  of  an 
incremental  change  of  the  output  voltage  with  respect  to 
variations  in  one  resistive  element  is 

AV  (t,R) 

— - -  =  I  (t)*I  (t)  (4.30) 

AR 


If  it  is  required  to  find  the  variation  of  V 2  with  respect 
to  all  resistive  parameters  the  changes  are  added  due  to  the 
superposition  principle  ever  all  branches,  a,  containing 
resistors.  Therefore 


AV  (t,R) 

— - -  =  Zi  (t)*i  (t)  (4.31) 

AR  a  a  a 

For  inductive  parameters  =  sL  and  x  =  L  AV2(t,L)  is 
determined  by  the  time  derivative  of  the  convolution  between 
the  corresponding  inductor  currents  in  the  original  network 
and  its  adjoint.  Therefore 


AV2(t,L) 

AL 


yd_,  ■ 
adt  Ua 


(t) *iQ  (t) ) 


(4.32) 


Finally  for  capacitive  elements  the  sensitivity  of  V2 
due  to  changes  in  all  capacitors  turns  out  to  be  the  summation 
over  all  capacitive  branches  of  the  time  derivative  of  the 
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convolution  between  the  voltages  across  corresponding  capacitors 
in  the  original  and  its  adjoint  network.  Therefore 


=  zlr  (V  (c)*V  (t)) 
AC  dt  a  a 


(4.33) 


a 


Starting  with  equation  (4.29)  the  computational  process 
can  be  simplified  by  considering  the  adjoint  network  (Fig.  8) , 
excited  by  a  current  source,  I 2,  as  follows 


9u 

X  =  ( — — ) I  (4.34) 

z  3x  a 


The  current  through  the  variable  parameter,  I  ,  is  then  given 
by 


T 

"a 


9V2(s,x) 

9x 


(4.35) 


Remembering  that  E  (s)  is  unity  in  the  frequency  domain  and  a 

9 

unit  impulse  in  the  time  domain,  the  sensitivity  of  the  output 
voltage,  AV2,  is  given  by 


AV~(s,x) 

— - -  =  I  (s,x)  (4.36) 

Ax 

in  the  frequency  domain,  and  by 


AV„  (t  ,x) 

— - -  =  i  (t,x)  (4.37) 

Ax  a 


in  the  time  domain. 

If  the  interreciprocity  theorem  is  applied  to  the  circuit 
of  Fig.  8,  interchanging  excitation  source,  I2,  and  response, 
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Iq/  the  sensitivity  model  is  obtained  as  shown  in  Fig.  9a. 

Fig.  9a  holds  for  the  frequency  domain  as  well  (t  replaced  by 
s) .  This  derivation  is  valid  for  all  types  of  passive  network 
parameters . 

For  the  transition  of  the  dependent  sources  the  voltage 
dependent  voltage  source  is  chosen ,  where 

-  ua(s;x)  =  Voy(s,y)  (4.38) 


such  that 
x  =  y 


Equation  (4.29)  becomes 


AV2(s,y) 


V.  I 
iy  oy 

E“TsT“ 

g 


5Voy(s'^ 

3y 


Ay 


(4.39) 


From  Fig.  8  the  excitation,  I2,  becomes  a  voltage  source  of 
value 


*2  = 


aVou(s^ 

3M 


(”V 


(4.40) 


Then  the  current  in  the  dependent  branch  of  the  current  dependent 
current  source  in  the  adjoint  network,  I  ,  gives  the  desired 
sensitivity 


I  (s)  = 
oy 


AV2 (s,y) 
Ay 


in  the  frequency  domain  and 


i.  (t)  = 
iy 


AV2  (t , y ) 
Ay 


in  the  time  domain. 


(4.41) 


(4.42) 
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Figure  9.  Sensitivity  Models  Derived  by  Application 
of  Tellegen's  Theorem  and  the  Mutual 
Recprical  Adjoint  Network 

a)  For  One-Port  Passive  Parameters 

b)  Voltage  Dependent  Voltage  Source 
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Interchanging  the  excitation  and  the  response,  referring 
to  Fig.  8  and  Fig.  9a,  leads  to  the  sensitivity  model  of  a 
voltage  dependent  voltage  source.  The  excitation  voltage  source, 

*v* 

given  in  (4.40)  has  to  be  transferred  into  the  time  domain 
such  that 


ig (t)  “ 


- *(-V.  (t)) 


(4.43) 


Knowing  that 


y 


+ 


vi 


(4.44) 


gives 


i2(t) 


("3ir  +  V‘<-Viu(t)) 


(4.45) 


The  minus  sign  in  front  of  V..  means  that  the  excitation  source 

l  M 

in  the  sensitivity  model  has  to  oppose  the  normal  current  flow. 
To  be  consistent  with  the  structure  of  a  dependent  source  and 
the  equation  (4.45),  V.  has  to  be  as  follows 


(t) 


9v.(t) 

9y 


(4.46) 


Substituting  (4.46)  into  (4.45)  leads  to  the  sensitivity  model 
(Fig.  9b)  as  stated  in  Chapter  III. 


C.  COMPARISON  OF  THE  TWO  APPROACHES 

As  a  main  conclusion  it  can  be  stated  that  the  sensitivity 
model  is  not  an  independent  method  for  computation  of  network 
sensitivity  but  a  special  case  of  the  mutual  adjoint  network. 
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For  sensitivity  calculations  in  the  frequency  domain  with 
respect  to  single  element  perturbations  both  methods  are 
equally  well  suited.  If  the  total  increment  in  sensitivity 
due  to  variations  in  several  parameters  is  required,  the  adjoint 
network  approach  is  advantageous  because  still  only  a  single 
excitation  is  required  and  the  analyses  of  only  two  networks 
at  each  frequency  point  are  necessary.  This  is  in  contrast 
with  the  use  of  the  sensitivity  model  where  a  separate  source 
is  required  for  each  variable  element.  This  requires  the 
analysis  of  one  network  for  each  parameter  at  each  frequency 
point. 

In  the  time  domain  the  sensitivity  model  is  the  better 
approach,  especially  if  single  parameter  changes  are  involved. 
The  desired  sensitivity  requires  the  analysis  of  one  network 
only  and  the  answer  comes  out  immediately  in  the  time  domain. 

In  contrast,  the  adjoint  network  approach  involves  the  analysis 
of  two  networks  and  requires  convolution  of  the  corresponding 
circuit  responses.  Alternately,  in  the  time  domain,  the  adjoint 
network  may  be  excited  by  a  source  with  backward  running  time. 
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